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Integrable flows on the Grassmannians Gr(N-l,N+l) are denned 
by the requirement of closedness of the differential N-l forms fijv_i 
of rank N-l naturally associated with Gr(N-l, N+l). Gauge-invariant 
parts of these flows, given by the systems of the N-l quasi-linear differ- 
ential equations, describe coisotropic deformations of (N-l)-dimensional 
linear subspaces. For the class of solutions which are Laurent poly- 
nomials in one variable these systems coincide with N-dimensional 
integrable systems such as Liouville equation (N=2), dispersionless 
Kadomtsev-Petviashvili equation (N=3), dispersionless Toda equation 
(N=3), Plebanski second heavenly equation (N=4) and others. Gauge 
invariant part of the forms fijv-i provides us with the compact form 
of the corresponding hierarchies. Dual quasi-linear systems associ- 
ated with the projectively dual Grassmannians Gr(2, N+l) are defined 
via the requirement of the closedness of the dual forms f2jv-i- ^ ^ s 
shown that at N=3 the self-dual quasi-linear system, which is associ- 
ated with the harmonic (closed and co-closed) form Q2, coincides with 
the Maxwell equations for orthogonal electric and magnetic fields. 
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1 Introduction 



Multidimensional integrable systems which are equivalent to the commuta- 
tivity of multidimensional vector fields 

N d 
D a = ^2a ak (x; A)- — , a = 1,2 (1) 

k=0 k 

where a a k(x; A) are rational functions of the spectral parameter A have been 
invented by Zakharov and Shabat in the seminal paper [I]. Concrete ex- 
amples of such systems and an extension of this scheme to the case when 
operators D a contain also the derivatives with respect to spectral parameter 
A (Jx = ■£-) have been considered in the papers [2], pg], [51 El El, [HIE], 
[TUl [TT] . [12] . This class of multidimensional integrable systems includes 
second heavenly equation [2], hyper-Kahler hierarchy [31 0], dispersionless 
Kadomtsev-Petviashvili (dKP) equation El H]> Manakov-Santini system 
[SI [H], Boyer-Finley (dispersionless Toda) equation [TUl E], Dunaski equa- 
tion [12] and their generalisations. In recent years various properties of such 
integrable systems (wave breaking, reductions, etc.) have been intensively 
studied [Ml [HI HH HH HH [201 IZI|. In particular, in the paper [TJ it 

was observed that in addition to the equations for the variables a a k{x) the 
equations for the Jacobian J of the solutions \1/ of the systems D a ty = 0, 
a = 1, 2 play an important role in the whole method. For instance, it allows 
to get a compact form of the corresponding hierarchies. 

In the present paper we show that all these integrable systems arise nat- 
urally in the families of the Grassmannians Gr(N-l,N+l) parameterized by 
N+l variables Xq, x±, . . . , x^. Namely, they are equivalent to the closedness 
conditions for differential N-l forms 

Vt N _i = uj A • • • A lo N -2 (2) 

where 1-forms up = J^fJi 1 Pi^Xi, P = 0, . . . , N — 2 are built of N-l vectors p@ 
defining (N-l)-dimensional linear subspaces. In terms of Pliicker coordinates 
n i h ...ijv_ 2 ; ik — 0, 1, . . • , N one has the system of linear equations 



9TTi i± ...i N -2 

dxi N _ 1 



(3) 
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This system is equivalent to the system of N+l quasilinear equations 

dJ + y^ 2 d(Ja lm ) _ Q dJ + y^ 2 d{Ja 2m ) _ Q ^ 
dx N _i dx m ' dx N dx m 

daik da 2k da lk da 2k \ , AT , . 

a a + 2^ a ^ a u~a~ I = °> k = 0,1, N — 2. (5) 

for J = 7To ) i ) ... ) jv-2 an d 2(N-1) independent affine coordinates 

O'lk = ( — l) fc ^ l7r 0...fc-lfc+1...7V-l) a 2k = (~ 1)^^ l7r 0...fc-lfc+l...Af, 

where k = 0, . . . , N — 2 The subsystem (jSj) coincides with the system describ- 
ing the coisotropic deformations of the N-l dimensional linear space defined 
by the equations 

N-2 N-2 

Pn-i + ^2 aikPk = °' Piv + ^ fl 2fcPfc = 0. (6) 

fc=0 k=0 

The system fl3]) provides us with the linearization of nonlinear equation (JSJ) 
in terms of the variables J, Ja\ m , Jci2m, where J is a solution of equations 

©■ 

For the class of solutions polynomial in the variable xq(= A) the system 
(JSJ) becomes the system of the integrable N-dimensional equations of the 
type considered in the papers [2] , [31 g] , El [7] , [9] , [101 E] , [12], plj|20]. 
Hierarchies of systems arising in such a way can be represented in the compact 
form 

= (J" 1 ^ A A • • • A d^ N _ 2 )_ = 

where ^o, . . . , ^n-2 ar e Darboux type cooordinates for the closed N-l 
form Qn-i which are the Laurent series in xq = A and (•■■)- denotes the 
projection on the negative part of Laurent series. Simple generating formulae 
for the variables ty k are also presented in [20] . 

Duality between the Grassmannians Gr(N-l,N+l) and Gr(2,N+l) sug- 
gests to consider the system dual to the system fH|5|) . It is equivalent to the 
closedness condition for the 2-form = *Qn-i where * denotes the Hodge 
star (or duality) operation or to the co-closedness of the form Qn-i an d it is 
of the form 

|^ + E^ kl = 0, k = 2,...,N, (7) 
dxk ^ dx m 
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= 0, Z = 0,1; j,k = 2,...,N 



(8) 

where J* = n^, a* = J* 1 ix{ j ,a* l = -J* 1 7r^ and 7rL are Pliicker coordi- 
nates dual to 7Tj 0j j 1 ...i^_ a • At N=3 the system (jSJ) is quite similar to the 
original system (J5J). For N> 4 the system OH]) is the semi-decoupled system 
of ( JV ~ 2 K Af ~ 1 ) four- dimensional subsystems. 

Self-dual systems are associated with the harmonic forms fijv-i, i-e. the 
forms which are closed and co-closed. Self-dual system in Gr(N-l,N+l) is the 
systems of (N — l) 2 equations for 2(N-1) dependent variables aik,a 2 k,k = 
0,1,..., AT — 2. At N=3 such self-dual system coincides in form with the 
system of sourceless Maxwell equations with vanishing second invariant (or- 
thogonal electric and maghetic fields) that provides us with a wide class of 
their solutions. 

The paper is organized as follows. Basic facts on the Grassmannians are 
collected in section 2. Closed differential N-l forms and associated systems 
of quasi-linear differential equations are considered in section 3. Solutions of 
these systems which are Laurent polynomials in one variable and their rela- 
tion to N-dimensional integrable equations are studied in section 4. Compact 
form of these equations and corresponding hierarchies are discussed in section 
5. Dual and self-dual quasilinear systems are considered in section 6. 



2 Grassmannians Gr(N-l,N+l). 

Grassmannian manifold Gr(m,n) is by definition the parameter space for the 
totality of m-dimensional linear subspaces V m in the n-dimensional space V n 
(see e.g. [22], [22], [21]) • Equivalently 

Gr(m,n) = {m-frames in V n } /GL(m) (9) 

where an m-frame means an m-tuple {p ,^ 1 , . . . ,p m_1 } of linearly indepen- 
dent vectors with the coordinates (i=0,l,. . . ,n-l) in a given basis {e°, e 1 , . . . , e n_1 } 
in V n . These coordinates can be arranged in the n x m matrix P with the 
elements . The dimension of Gr(m,n) is equal to m(n-m). 

Grassmannain Gr(m,n) can be viewed as an algebraic submanifold of the 
( n _^) !m! -dimensional space A m V n via the correspondence of a m-frame {p°, 
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p 1 , . . . , p m 1 } with the exterior product p° A p l A ■ ■ ■ A p m 1 G A m V„ (the 
canonical projective embedding). One has 

p° Ap 1 A---AP™- 1 = ^oi 1 ... !m -/°Ae' 1 A---Ae t "- 1 (10) 

0<io<— <im-l 

where vr ioil ... im _ 1 = det(p|) fe) ; =0 ,..., m -i. These coefficients 7r io fl ... im _ 1 are called 
the Pliicker coordinates and satisfy the Pliicker's relations 

m 

E(- 1 )^o...^- 2 ^ ..jj...c = ° ( n ) 

where indices it and z'^ range over all possible values and the notation i\ 
means the deletion of this number. The Pliicker coordinates completely define 
embedding of Gr(m,n) into A m V n . 

Pliicker's coordinates define also the subspaces V m in V n . Namely, a 
point in V n with the coordinates (2/0,2/1, • • • , Un-i) lies in V m iff (see e.g. [22]). 
they obey the system of equations 

m 

^(-l)^ o ... i; „ lii+1 ... lm 2/,=0. (12) 
1=0 

In virtue of the Pliicker's relations there are n-m linearly independent equa- 
tions among them. As a consequence, the system ffTTj) is equivalent to the 
system of n-m equations 

m— 1 

yy + ^2 a ykVk = 0, 7 = rn,m + 1, . . . ,n - 1. (13) 

k=0 

where independent affine coordinates in Gr(m,n). 

The Grassmannian G(m,n) is projectively equivalent to the Grassmannian 
Grfn-m.n') 1221. The latter has coordinates it* ,■ defined by the relation 

do ii ...i N K m i m+1 ...i N = n io h ... lm _, , where e k) h ... ijv is the totally antisymmetric 
tensor in N+l dimensions. The coordinates n* m i m+1 i N obey equations ( TTTj) 
and the system dual to (fl2"]l is 

^<...^ mJ ^r = o (14) 
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This system of equations is equivalent to the following one 



n—m—l 



y*i + a *~f kV * k = °' l = n-m,...,n (15) 



k=0 



where a* k are independent affine coordinates in Gr(n-m,n). 

In this paper we will consider Grassmannians Gr(N-l,N+l) for arbitrary 
N. They are the 2(N-l)-dimensional sets of N-l- dimensional linear subspaces 
in the N+l- dimensional space defined by codimension two constraint 



7V-2 N-2 

y N - 1 + ^2a lk y k = 0, + a^kVk = (16) 

k=0 k=0 

with a lk = (-l) k J- 1 7i ...k-ik+i...N-i, a 2 k = (-l) k J~ 1 no...k-ik+i...N, k = 
0, . . . , N — 2, J = 7Toi...at-2- Affine coordinates a\ k and a 2k are completely 
defined by the set of N-l independent vectors in Vjv+i- 

In particular, at N=3 for the Grassmannian Gr(2,4) equations f[T5|) are of 
the form 

Vi + any! + a w y = 0, y 3 + a 2 \y\ + a 20 y = 0. (17) 
Pliicker's coordinates obey the single equation 

7T017T23 - 7T027T13 + TTq^TTu = (18) 

and 



(19) 

where p° = (Po, P®, P 2 j P3) an d P 1 = (Po, P11 P 2 , P3) are two vectors in V4 defin- 
ing the plane (fT6|) . Due to the relations (TIE]) one has 



vi"02 = — Jan, 7T 3 = - Ja 2 i, n\ 2 = Ja 10 , 7Ti 3 = Ja 20 , vr 23 = J(a 2 iai - 011020)- 

(20) 
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The Grassmannian dual to Gr (N-1,N+1) is Gr(2, N+l). It has the same 
dimension 2(N-1) as the former one , but is the set of 2-dimensional linear 
subspaces in the N+l-dimensional space Vat+i defined by the constraints 

y* + a* iy l + a; o y* = 0, j — 2, . . . , N. (21) 

At N=3 equations (ITT)) . (ITS)) and their dual equations ( 121]) have the same 
form. 

3 Closed differential N-l-forms and systems 
of quasi-linear differential equations 

Now we introduce an infinite family of Grassmannians Gr(N-l,N+l) parametrized 
by N+l variables x , x±, . . . , Xn- Thus, all quantities considered in the pre- 
vious section become functions of these variables. A geometrical realization 
of such family is provided, for example, by the 2(N+l)-dimensional cotangent 
bundle with the local coordinates xo, x±, . . . , x^- in the base manifold and 
cooordinates p , pi, . . . , p^ in cotangent space T^ +1 and the Grassmannian 
Gr(N-l,N+l) in attached to each point of the base space. 

For each vector p 13 defining the linear N-l dimensional subspace in V)v+i 
there is a naturally associated differential 1-form up = Y2iLoPi dx i ■ ^ ne 
projective embedding of the Grassmannian Gr(N-l,N+l) ( formula f fTOj) ) 
suggests to introduce differential N- 1-form 

^at-i = coo A • • • A (22) 

associated with the member of the family of Grassmannians corresponding 
to the point (x , Xi, . . . , xn). One has 

ftjv-i = ^2 *i - ^-2 ( x ) rfx *o A dx ii A ' ' ' A dx iN~2 ( 23 ) 

0<io<---<«JV-2 

where all indices take values 0,1,.. .,N. The homogeneous coordinates ■K iQil ... iN _ 2 
are defined up to the common factor p , depending on x. So, the form f )23|) is 
defined up to the scaling transformations f27v_i — > p(x)Qn-i with arbitrary 
function p(x) which can be viewed as the gauge transformations. 

To select a special family of Grassmannians Gr(N-l,N+l) we require that 
this form is closed. Thus, we require that the Pliicker coordinates obey the 
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equations 

dTTj j± ...ijv_ 2 _ q (24) 

dxi N _ 1 

where indices z& take values 0, 1, . . . , N and the bracket [ . . . ] means antisym- 
metrization over all these indices. This system of 2 J^" 1 = N + 1 differential 

°JV-1 

equations toqether with the algebraic Pliicker relations (ITT]) form a full sys- 
tem of equations which characterizes family of Grassmannians Gr(N-l,N+l) 
for which the form Ojv-i is closed. 

The system (j2"31) is equivalent to the system 



dJ + ^ 2 g( j fllm ) = + y g^gjnO = Q (25) 



<^7V-1 ^ & m ' dx N ^ dx 



da lk da 2k ,V^V daik da2k \ n i. n i jvr o fo^ 

+ ~ ail ~dx~i) = ' = 0, 1, N — 2. (26) 



The subsystem (1261) . invariant under the gauge transformations mentioned 
above, is the gauge invariant form of the system ( l24j) . Two equations ( 125]) 
can be viewed as the equations for the gauge variable J which transforms 
as J — > pJ under the gauge transformations. It is a straightforward check 
that the equations (l25l) are compatible due to the subsystem (126]) . We note 
that the system ( f25|26]) can be rewritten in the form 

N-2 o N-2 „ 

O.W + i:^= = 0, fl2 lnJ + ^|^ = 0, (27) 

m=0 (7Xm m=0 (7Xm 

£> 2 a lfe - Dia 2fc = 0, fc = 0,...,JV-2 (28) 
where -Di and D 2 are the vector fields 

^ O \ ^ C _ C \ O /_„\ 

+ >aimo— , ^2 = tt— + > a 2 m.T^. 29 

OTiV-l ^ OXN ^ 

So the subsystem ( 126]) is equivalent to the commutativity [Z?i, -D 2 ] = of the 
vector fields D\ and D 2 . 

The system (126]) of nonlinear equations admits an obvious linearization. 
Indeed, for a solution a±k a 2 k of this system one can find J solving linear 
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equations (1251) Then for variables J, Ja^Ja^ one has a linear system 
which coincides with the system (I24p . 

The system ( 126]) has another geometrical meaning Let us denote by J the 
family of ideals of the N-l-dimensional linear spaces defined by the equations 
ffTB"]) and require that J is the Poisson ideal, i.e. 



{J, J} C J 



(30) 



where { , } is the standard Poisson bracket with Xo, Xi, ■ ■ ■ , x^] Ho, Hi, ■ ■ ■ ,Un 
being the canonical Darboux coordinates. Condition f[3"0"j) defines the so- 
called coisotropic deformations of the linear algebraic variety defined by 
equations (jl~2]) [25j,[26j. Such deformations are described exactly by equa- 
tions (1261) . Thus, the coisotropic deformations of the ideal J represent the 
necessary and sufficient gauge- invariant conditions for the closedness of the 
(N-l)-forms (122]) . 

At N=2 we have two-dimensional family of straight lines in the 3-dimensional 



space, 1-form f2i 
by equations (J 



^2 k=0 Pkdxk and the conditions of its closedness are given 
Po) 



dp dpi 



dx\ dxc 



0. 



dp dp 2 



dxr> dxv, 



0. 



dpi dp 2 



dx 2 dxi 







or 



dJ | d(Ja 
dxi 
da w 



10, 



dx 
da 20 



0. 

da 



dJ d(Ja 20 ) 



dxo 



dx-2 



dxi 



«20 
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dx 
J- 1 - 



aw 



dx 
da 2Q _ 
dxn 



0. 



0. 



(31) 



where J = poa 10 = - J 1 p 1 , a 20 = -J ^p 2 . 

At N=3 one has family of planes (TXT]) . 2-form Q 2 given by 

fl 2 = J(dxo A dx\ — audxo A dx 2 — a 2 \dxo A dx^ + a^dxi A dx 2 +(32) 
+a 2 odxi A dx3 — (aua 20 — ai a 2 i)dx 2 A dx^) 



and the equations 
dJ 



_ + d{Ja lr 
dx 2 ^ 

m=0 



dx r 



0. 



dJ ^ 

dx 3 ^ 



m=0 



d(Ja 2r , 



(33) 
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da lk da 2k , / da lk da 2k \ n , „ . /0/A 
£(^^-«i^J=0, fc = 0,l- (34) 



dx 3 dx 2 



Finally, we note that at arbitrary N the N-l form Qn-i has a factorized 
form 

= J&N-1 (35) 

where is the gauge invariant N-l-form depending on the affine coordi- 

nates etifc, a 2k only. 

4 Solutions rational in one variable and N- 
dimensional integrable systems 

The system ( 12"6"|) of N-l equations for 2(N-1) dependent variables always is the 
underdetermined system . In order to make it determined one should imposes 
constraints. The simplest constraint is to choose a 2k to be certain functions 
of ai fc , i.e. a 2k = Pk(aio, an, . . . , ai N - 2 ), k = 0, 1, . . . , N — 2. In this case 
the system (26) becomes the system of N+l-dimensional hydrodynamic type 
equations 

^ + E^F^ + E^V^ = ' * = 0,l,...,tf-2 (36) 
ox N ^ dx N -! ^ dxi 

where 

A km = —jt^-, B kml = 5 km ipi - au-^^, k, m, I = 0, 1, . . . , N - 2. (37) 
At N=2 it is effectively the 1+1- dimensional Hopf-Burgers type equation 

~dT ^ ~dxti ~ (38) 

where A—l — J^, B = ip - aio J£ and t = \{x 2 - x{). 

Different types of reductions of the system (1261) are associated with the 
restriction to a special subclass of its solutions, for example, to solutions 
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which are Laurent polynomials in one variable, say Xq. Considering, for 
example, solutions of the systems ( 12 6 p of the form 



a lk = ^ U kn(Xl, X N )Xq, a 2k = ^2 v kn{x\, • • • , X N )Xq, (39) 



n=—m\ n=—m2 



where k = 0, 1, . . . , N — 2, with appropriate ni,mi,n 2 ,m 2 , one gets the 
systems of N-dimensional equations for the functions v kn together with 
the corresponding constraints between them. A complete classification of 
determined systems which can be obtained in such a way is beyond the 
scope of this paper. Here we will present several examples to demonstarte 
that one can get number of known integrable systems within this approach. 
In the simplest case N=2 with the Ansatz 

oio = u + Aui, a 20 = v + Ai>i + X 2 v 2 (40) 

where A = x , x = (xi, x 2 , x 3 ) the system f[3"Tj) assumes the form 

u 0xa -v OsB1 +u 1 v -uoV 1 = i u lx2 -v lxi -2u v 2 = 0, v 2xi +u 1 v 2 = 0. {±1) 

This system implies that u\ = —(p xi where <p = \nv 2 . Under the constraint 
Vi = the above system becomes 

<p xixa + 2u Q e v = 0, u 0x2 - v xi e-* = (42) 

where v = 2v e ip . Under the further constraint u Q = —\,v = 2 one gets the 
Liouville equation 

(Pans* = (43) 



for which 



Choosing 



aio = ~\ ~ A^ 1; a 20 = + AV. (44) 



a w = -- - Xip xi , a 20 = e 9 + X 2 ip xlX2 + A 3 e 2</ \ 



2 

one obtains the higher Liouville equation 



3 

¥x lXl x 2 - Vx 1 ^x 1 x 2 - = ( 45 ) 
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which in terms of the variable \I/ = exp(— ^<p) is of the form 

_ ^-4 = (46) 

Solutions of the Liouville and higher Liouville equations describe two- 
parameter deformations of the straight line defined by the linear system 

Pi + aioPo = 0, p 2 + a 20 p = (47) 

in the three-dimensional space with a\oa 2 Q given above. In geometrical 
terms solutions of these equations describe special classes of congruences of 
the lines in the three-dimensional space. In a different context the interrela- 
tion between congruences of lines and system of equations of hydrodynamic 
type has been studied in [27]. [28]. 

In the case N=3 the first nontrivial choice is 

a 10 = u (x), au = u 1 (x)+\, 

a 2 o = v (x) + Xvi(x), a 21 = v 2 (x) + Xv 3 (x) + X 2 . (48) 

Substituting this ansatz into two equations (|34|) and denoting x = xi, y = x 2 , 
t = x 3 , one gets the system of equations 

u ot - v 0y + v 2 u 0x - uiv 0x - uovx = 0, (49) 

U\ t - V 2y + Vq + V 2 U lx - U X V 2x - UqUi = 0, (50) 

vi - 2u - v 2x - v 3y + v 3 u lx - Uiv 3x = 0, (51) 

vi y + v 0x + UiV lx - v 3 u ox = 0, (52) 

uix ~ v 3x = 0, (53) 

u 0x ~ v lx = 0. (54) 

Equations (!53l) . (|54l) imply that V\ = Uo t v 3 = u\ and equations (loTI) . (1521) 
become 

uq + ui y + v 2x = 0, u 0y + v 0x = 0. (55) 

Last equation (15151) implies the existence of the function u such that n = 
— u x , Vq = u y . As a result from the first equation (151)1) one gets Ui = v x , v 2 = 
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u—Vy, where v is some function. Substituting these expressions into equations 
(1501) . one obtaines the system 



U xt + Uyy +Ul+ (U — Vy)U XX + V X U X y = 0, (56) 

v x t + v yy + v x v xy + (u - v y )v xx = 0. (57) 

It is the Manakov-Santini system found in [8] , [9J . Under the constraint v — 
it reduces into the dKP equation 

u x t + Uyy + [uu x ) x = (58) 

while at u = it is the equation 

V X t + Vyy + V x V xy - VyV xx = (59) 

considered in [29j [30], [31]. Thus, the Manakov-Santini system is the gauge- 
invariant form of the closedness condition for the 2-form 

Q 2 — J(d\ A dx — (v x + X)dX A dy — (u — v y + \v x + \ 2 )d\ Adt — u x dx A dy+ 

+ {u y — \u x )dx Adt — [u x u — u x v y + u y v x + Xu y )dy A dt) (60) 

associated with the family of Grassmannians Gr(2,4). The function J is a 
solution of the system 

J y ~ u x J x + (A + v x )J x + v xx J = 0, (61) 

J t + (u y - Xu x )J x + (u — v y + Xv x + A 2 ) J x + {-v xy + Xv xx )J = 0. (62) 

Second example of the 3-dimensional integrable system corresponds to 
the choice (again A = xq, x = x\, y = x 2 , t = £3) 

aio = -K4>y - 4>x— ), an = -A - — , 
m x m x 

.<f> x 1 e* , s 

fl2o = -e. — , a 2 \ = —r — , (63) 
m T X m T 



where (ft and m are functions of x, y, t. The system (1M|) takes the form 

{e*)xx ~ fn x (p t y + m y (f) tx = 0, (64) 
e*Tn xx - m x m ty + m y m tx = 0. (65) 
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It is the generalization of the 2DTL given in [2T] . Under the reduction m = x 
it is the 2DTL equation 

(e% x -<f>t v = (66) 
and at = it is the equation 



m. 



- m x m ty + m y m tx = (67) 

considered in [29], [31], [30]. 

The corresponding 2-form is given by 

Tfl 1 6^ Tfl 

Q 2 = J(dX A dx + (A -)dX Ady+ - — dX A dt + A(0„ - <p x —)dx A dy+ 

m x X m x m x 

+ e^—dxAdt (A0 X + d> v )dy A dt) (68) 

m x m x 

and J is a solution of the system 

J y + X(<f> y - 0,^) J X - (A + —)J X + (0, - 0,^ - (^)„) J = 0, 

e*0 x 1 1 , e^ 1 N 

Jt - — </a - -r J x ~ T — «J = 0. 69 

Now let us consider the case N=4 and choose the functions a±k and a 2 k, 
k — 0, 1, 2, as 

a 10 = ii , an = ui + A, a 12 = u 2 , 

020=^0, «21 = Ui, a 2 2 = ^2 + A (70) 

where t^; = 0, 1, 2 are functions of the variables x = x±, y — x 2 , z = X3, 
t = X4 and A = xq. Substituting this ansatz into the system (1261) with N=4, 
one gets the equations 

Uot ~ VQz + V\Uq x + V 2 U 0y - U 1 V Qx - U 2 V 0y = 0, (71) 
Uu - Viz +V + V!U lx + V 2 Ui y - U!V lx - U 2 V ly = 0, (72) 

u 2t - v 2z - u + v x u 2x + v 2 u 2y - u x v 2x - u 2 v 2y = 0, (73) 
Uk y -Vkx = 0, A; = 0,1,2. (74) 
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The last equations imply that 

U = U x , V = Uy, Ul = V x , Vl = Vy, U 2 = W x , V 2 = Wy (75) 

where w, u, v are functions obeying the equations 

Uxt - Uyz + v y u xx + (w y - v x )u xy - w x u yy = 0, (76) 

Vxt - Vy Z + Uy + VyV XX + (Wy ~ V X y ~ W X V yy = 0, (77) 
Wxt - Wy Z ~ U x + VyW XX + (Wy ~ V X )W X y ~ W X W yy = 0. (78) 

The 3-form Q3 is the sum of 10 terms 

Q3 = J(d\AdxAdy—w x d\AdxAdz—(\+Wy)d\AdxAdt+(\+v x )d\AdyAdz+. . . ) 

(79) 

and equations for J are 

J z + u x J\ + (A + v x )J x + w x J y + (v x + Wy) X = 0, (80) 

J t + UyJ X + VyJ x + (A + Wy) Jy + (V x + Wy)y = 0. (81) 

The system (I76H78P admits several reductions. An obvious one is u = 
for which it is the system of two equations (1771 1781) with u = 0. Less trivial 
reduction corresponds to the constraint v x + w y — 0. It implies that v = Q y 
and w = —Q x . and the system (17611781) takes the form 

U X t Uy Z ~t~ QyyU xx 2@ X yU X y ~\~ O XX Uyy 0, (^^) 

u = Q zy - Q tx + Q 2 xy - Q xx Q yy . (83) 

It is the Dunajski system proposed in [12]. Under the additional constraint 
u = one has the Plebanski second heavenly equation [2]. 

<3> zy - <3>t x + Qly - <3> xx <3> yy = 0. (84) 

Integrable systems for any N can be constructed in a similar way. An 
interesting problem of finding the systems which describe solutions of the 
system (I26p which have poles in one variable, similar to the Calogero-Moser 
system for the rational solutions of the Korteweg - de Vries equation [32] , 
will be discussed in a separate paper. 
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5 Compact form of the hierarchies 



Considering the functions a\k and a 2 k which are higher order polynomials in 
xq = A, one gets higher Manakov-Santini and Dynajski equations (for sim- 
plicity we will not consider here two-component 2DTL case which requires 
Laurent polynomials, however, it can be also considered in the similar frame- 
work, see [211 EE])- The forms fijv-i provide us also with a compact forms 
of these hierarchies. These differential forms have rank N-l. Hence, the con- 
dition cIQn-i = implies the existence of N — 1 variables \l/ , ty 1 , ty N ~ 2 
such that (see e.g. [33| IM]) 

Q N _ 1 = A d^ 1 A . . . cM N - 2 . (85) 

This means that the components of the vectors p k , k = 0, 1, . . . , iV — 2 which 
define N-l- dimensional linear subspaces in the Grassmannian Gr(N-l,N+l) 
can be taken as the derivatives p k m = k — 0, 1, . . . , N — 2; m — 0, 1, . . . , N 
and equations f fT6|) take the form 

£>i** = 0, D 2 ^ k = 0, fc = 0,l,...,JV-2 (86) 

where the operators D\ and D 2 are given by (T29]) . 

Since fi/v-i is a complicated function in xq = A the variables \l/ fc are, in 
general, certain Laurent series in A. An important property of this form is 
that the N-l-form Cl^_i defined in (155]) is a polynomial function in A for 
polynomial a\k and a 2 k- Hence, 

(j-WAtf A...^- 2 )_ =0 (87) 

where (• • •)_ denotes the projection on the part of (• • •) with negative powers 
in A and 

J = 7T l...jV-2 = det(<9^ m ) mi/=0 ,...,7v_2- 

This is the compact form of the hierarchies of multidimensional integrable 
systems considered in [20] . 

Generating relation ( 187]) implies Lax-Sato equations which define the evo- 
lution of the series W 1 , . . . , ty N ~ 2 with the coefficients depending on x m (1 < 
m < N — 2) with respect to the times t\ = xat-i, h = %n (and also higher 
times corresponding to commuting flows of the hierarchy). Indeed, the coef- 
ficients of the form ( 155]) can be written as n io ix ...i N _ 2 = det(<9 i; \l/ m ) mi / = o i ...,Ar_2 
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(compare with relation f llOp ). The functions satisfy identically the rela- 
tions corresponding to formula (TT21) . 

N-l 

E(- 1 )^o...^ lii+1 ... JJ v^^ = 0. (88) 
z=o 

An important step is to observe that due to generating relation (l87|l the 
affine coefficients of linear relations (!88|) are analytic, so we obtain nontrivial 
relations 

m 

^(-l)V~V.i,- 1 * !+1 -J+^ = 0. (89) 

1=0 

Taking the basic relations corresponding to constraints f)16p . we obtain Lax- 
Sato equations 

N-2 

d N ^ + J2 a ikd k ^ q = 0, (90) 

fc=0 
N-2 

d N ^ + J2^kd k ^ q = (91) 

fc=0 

with 

Ctlk = ( — ^) k (J 1?r ...k-lk+l...N-2N-l) + , 
a 2k = ( — l) fc (^ l7r 0...fc-lfc+l...Af-2A r ) + , 

where k — 0, . . . , N — 2. Under reasonable assumptions the coefficients a±k 
and a,2k contain only finite number of the coefficients of the series ty q , and 
the equations ( )90|) . (|9Tj) define the evolution of these series with respect to 
times ti = Xn-i, £2 = xn (similar equations can be written for the hierarchy 
of commuting flows). The compatibility conditions for equations ( 190|) . (I9TT) 
considered as linear equations for \& 9 are of the form (1261) . and they define 
N-dimensional integrable systems for the coefficients of polynomials a^, a,2k- 
They can be also obtained directly from the condition of the closedness of 
the form Q. 

The general multidimensional integrable hierarchy arising from generating 
relation of the form flHTl) was considered in [20]. Here we give a brief out- 
line of this construction, emphasizing the connections with the construction 
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developed in the present work. First, following [2D], we impose generating 
relation (j87p on the formal series of the variable A 



v = \ + Y,K(t l , 

n=l 

oo 

** = j^t*(*°r+X;« r J;(t 1 ) ... ) 0(* 



0\-n 



(92) 
(93) 



n=0 



n=l 



where 1 ^ k ^ N — 2, depending on k infinite sequences of independent vari- 
ables t k = (£q, . . . , t*, . . . ), = Xk, X — Xo. Some motivation of introduction 
of the series of this type and higher times of the hierarchy is given in [3] 
in the simpler case of hyper-Kahler hierarchies (vector fields don't contain 
a derivative over the spectral variable). The setting of the hierarchy with 
infinite number of higher times corresponds to Gr(iV — 1, oo). 

Slightly modifying notations in the formulae (19"U|) . (191 p . we obtain infinite 
set of Lax-Sato equations of the hierarchy in the form 



N-2 



(94) 



k=0 

where I = 1, . . . , N - 2, n = 1, . . . , oo, & n = 



Jk 



■ir j 



d(^°,^ 1 



N-l -qjN-2 



D(x , ...,x k ,.. .,x N - 2 ,t l n ) 



and the fraction denotes a Jacobian matrix, x^ is the absent element. Using 
the series (19~2"|) . f[9"3"j) and estimating the projections, it is possible to simplify 
these expressions and get rid of the derivatives over higher times, 



-(-lf(-l) J~ (* 



-l/,TrO\n 



N-2\ 



D(x , . . . , Xk, • • • , XN-2) 



The corresponding Lax-Sato equations are 



d l n w = ^(_i) fc (_i)' j^^y 



k=0 



D(x , 



•Ek) 



%N-2 



These equations have an evolutionary form and define the dynamics of the 
series with the coefficients depending on x±, . . . 1X^-2 with respect to the 
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higher times. The evolutionary form of Lax-Sato equations is connected with 
the special choice of the form of the series (19~2"]1 . (19"H|) . it is a serious argument 
in favour of this choice. 

Using the Jacobian matrix 

(Jaco) = U(*o!...!*„_ 2 ) J ' det(Jaco) = J ' 

we obtain Lax-Sato equation in the form introduced in [20] . 

N-2 

dn* = E (( Jac or%(*T))+ W 1 ^ A; ^ iV - 2, (95) 

i=0 

where 1 ^ n < oo, * = (*°, . . . , * Ar " 2 ). It was proved in [20] that a complete 
set of Lax-Sato equations (|95|) is equivalent to the generating relation (157)1 
and that Lax-Sato flows are compatible. First flows of the hierarchy read 

Af-2 

= (Xd k - Ys( d * u p) d P - (d k uo)dx)V, (96) 
P =i 

where «o = w& = l^/c^A^ — 2. A compatibility condition for any 
pair of linear equations (e.g., with 9^ and 9', k ^ q) implies closed nonlinear 
N-dimensional system of PDEs for the set of functions u k , Uq, which can be 
written in the form 

d^d q u - dfd k u + [d k u, d q u] = (d k u )d q - (d q u )d k , 

d^dqUo - dfd k u + (d k u)d q u - (d q u)d k u = 0, (97) 

where u is a vector field, u = Yl p =i u k d k . For A^ = 4 this system corresponds 
to the system (1761175]) and, under a volume-preservation reduction J = 1, to 
the Dunajski system (152"] 155]) . The case A^ = 3 corresponds to the hierarchy 
connected with the system (156 ] l57j) (the Manakov-Santini hierarchy, see [20]). 

6 Dual and self-dual quasi-linear systems 

For the dual Grassmannian Gr(2,N+l) one defines the differential 2-form 

1^2 = ^n-i A oj n (98) 
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where 1 -forms o; 7 = J2f =0 pPdxi. So, 

^ = E<fl^oA^ (99) 

This 2-form and N-l-form Qn~i are connected by the Hodge star (duality) 
operation * 

tt* 2 = *£l N -i (100) 
where by definition (see e.g. [31]) 

*(dxi A dxi x A ■ ■ ■ A dxi m ) = 
= (n — m — 1 \\ ^ ] e io,!i,...,i m ,! m+ i,...,!„^ :j; ! m+ i A dxi m+2 A • ■ ■ A c?Xj n . (101) 

lib J. J ■ 

Operator 5 dual to the exterior differential g? is defined as 

5Q = {-l) mp+m+1 *d*Q (102) 

where m is the dimension of the space and p is an order of the form Q. 
In particular, 5d* = -kd8 and the operator A = d5 + 5d is self-dual. The 
differential form Q obeying the condition 5Q = ( or dQ* = ) is called 
co-closed. The form which is closed and co-closed, i.e. AQ = is refered as 
the harmonic form (see e.g. [S|). 

The condition of closedness of the form Q 2 1S given by the system of 
equations 



dxi 



(103) 

where indices take all values 0,1,.. .,N. In virtue of ( 12 ip this system is equiv- 
alent to 

dJ* d(J*a*) d(J*a*) 

oj + _A_Z^ + J_^lZ = , 7 = 2 ,...,iV, (104) 

ox 1 OXo OX i 

J m—U 

(105) 

At N=3 the system ( I104p . ( I105p consists from four equations which in the 
variables J* and a* t have the form similar to the original system (|33|) . (|34|) . 
In the original variables, due to the relation n* k = Ylim=o e ikim^im, i-e. 

tt^ = 2-7T23 = -2J(a 2 ia 10 -a 11 a 2 o), Hq 2 = -2n 13 = -2Ja 20 , Hq 3 = 2n 12 = 2Ja w , 
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n{ 2 = 2n 03 = -2Ja 2 i, 7T* 3 = -27T 02 = 2Ja u , ir 23 = 2n 01 = 2.1 (106) 
it is of the form 



dJ d(Ja 2 i) d(Ja n 



dx\ dx 3 8x2 
8J d(Ja 20 ) 8(Ja w ) 



0, (107) 
0, (108) 



8(Ja n ) 


d{Ja w ) 


d 


dxo 


dxi 


dx 3 


d(Ja 21 ) 


d(Ja 20 ) 


d 


dx 


dxi 


dx 2 



dxo dx 3 dx 2 

■[J(a 21 a w - a u a 20 )] = 0, (109) 

[J(a 21 a w - a n a 20 )] = 0. (110) 

Due to the equations (I107[ 11081) the equations (I109| II 101) are equivalent to 
the following 

da 10 da 1± da n da u da w da w 

a ^ fl io-^ ^ a 20^ an— a 21 — — = 0, (111) 

OXi OXq ox 2 OX3 ox 2 0x3 

da 20 da 21 da 21 da 2X da 20 da 20 

~a -r. ^«io^ ^ a 2o^ an— a 21 — — = 0. (112) 

OXi OXo OX 2 OX3 OX 2 OX3 



This dual system is quite similar to the original system 

If one requires that the form Q 2 is a harmonic one, then one has system of 
four equations (|34"|) . fillip , fll 12j) for four dependent variables aio, an, 020, «2i- 
This system has a simple physical meaning. Indeed, the conditions of closed- 
ness and co-closedness of the form D, 2 , i.e. the equations (J3J) and (11031) are 
equivalent to 

d7T ik dTTki dllu S^dTT ik 

—— + -— + -— = 0, > - — = 0, t, k, I = 0,1,2,3. (113) 
OXi OXi ox k ox k 

This system can be viewed as the second and first pairs of the source- 
less Maxwell equations for the electromagnetic field tensor 71^ , i.e. 7r 0a = 
E a , ir a/ 3 = — 5^ 7 =i e a/3-yH-y, a, (3 = 1, 2, 3 where e a/ j 7 is totally antisymmetric 
three-dimensional tensor. The Plucker's relation (18) in this terms is of the 
form 

3 

Y,E a H Q = 0. (114) 

a=l 

Thus, for the harmonic two-form Q 2 the equations for the Plucker's coor- 
dinates coincide with the sourceless Maxwell equations with perpendicular 
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electric and magnetic fields. This observation provides us with the variety 
of solutions of these equations (see e.g. [35]) Given such a solution one gets 
solution of the system fl3J}, ffTTTj) . ffTT2]) via 

H 3 E/2 H2 E 3 

£/l £/i £/i £/l 

and J = E\. 

It is a simple check that under the Manakov-Santini Ansatz (j48p the 
system (I34p . f illip . (I112p has only trivial solution. Rational solutions of this 
system will be studied elsewhere. 

At N=4 the dual system (11041) . (1 105[) is the system of 9 equations 



8J* d(J*a* ) 8{J*a* 2 



21 1 



8x2 dxo dx\ 
8J* , 8(J*a* 30 ) , 8(J*a* 31 ] 



8x3 dxo dx\ 
8J* d(J*al ) d(J*a* 41 ] 



8x4 8x dxi 



0, (116) 

0, (117) 
0, (118) 



da* 20 


8a* 30 


dx 3 


8x2 




8a* 31 


dx 3 


8x2 


da* 20 


8a* m 


dx^ 


8x2 


da* 21 


8a* Al 


8x4 


8x 2 


8a* 30 


8a* m 


8x4 


8x 3 


8a* 31 


8a* Al 



+ ^-^ + ^-^^0, (119) 

+ 030 8x Q ° 20 8x + %1 dx x 21 dx x ~ ' 1 ' 

+ «lo^ - + _ a * 2l ^k = 0, (121) 

8xq 8xq 8x1 8x\ 



MO 



<9xo cteo 8x\ 8x\ 



a »a^- a »a^ + a «a?- a ^ = - (123) 

, _ o . Saji + Soli _ . = „ (124) 
(9x4 ' cteo cteo ctei 9xi 

It is quite different from the original system ([231 |2"B"|) at N=4. The equa- 
tions (j!16| 11171 11191 1120p form a closed subsystem for the dependent vari- 
ables J*, aJjo, «2i5* a 30' a 3i- The same is valid for the groups of variables 
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J*, a2 , a^,* al , a* A1 and J* , a^ , a^* a* m , a* A1 . So, the whole system (11161 
I124p is decomposed into three independent subsystems. Each of these sub- 
systems coincides with the system ( |33]|3"4"|) for Gr(2,4). Such a decomposition 
remains valid also for the rational solutions discussed in section 4. So, the 
simplest rational solutions of the system ( I119til24p can be constructed as the 
common solutions of three independent Manakov-Santini or DToda systems. 
The subsystem (111911124]) . which is the gauge invariant form of the system 
( 11031) . describes the coisotropic deformations of the family of the planes de- 
fined by (l2Tj) . It contains twice more equations than the subsystem (126!) at 
N=4. 

Similar situation takes place for the dual systems for Gr (N-1,N+1) at 
N> 5. 
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